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We have introduced an improved X-ray phase-retrieval method with unprecedented speed of con- 
vergence and precision, and used it to determine with sub- Angstrom resolution the complete atomic 
structure of an ultrathin superconducting bilayer film, composed of Lai.55Sro.45Cu04 and La2Cu04 
neither of which is superconducting by itself. The results show that phase-retrieval diffraction tech- 
niques enable accurate measurement of structural modifications in near-surface layers, which may be 
critically important for elucidation of surface-sensitive experiments. Specifically we find that close 
to the sample surface the unit cell size remains constant while the copper-apical oxygen distance 
shows a dramatic increase, by as much as 0.45 A. The apical oxygen displacement is known to have 
a profound effect on the superconducting transition temperature. 
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The exciting discovery of interface superconductivity 
in complex oxides [TJ |2j |3j IH |5] has triggered intense de- 
bate about its origin and the possibility to enhance Tc 
even further [BJ |71 HI IH]- The interfacial enhancement 
[3] of the superconducting critical temperature (Tc) is 
influenced by the crystal structure. The Z-axis lattice 
constant (cq) varies significantly among La2Cu04 (LCO), 
Lai.55Sro.45Cu04 (LSCO), and bilayer LSCO/LCO films, 
depending even on the deposition sequence, and it affects 
superconductivity: Tc scales with Cq almost perfectly 
linearly [TD]. The reason for this is not understood at 
present, but notice that in (La,Sr)2Cu04 the change in 
Co goes together with the change in ca, the distance be- 
tween copper and the nearest apical oxygen, which some 
believe to play a key role in the high temperature super- 
conductivity (HTS) phenomenon dH [H [H [H [13 US] . 
In any case, it is certain that (i) from one cuprate to an- 
other, Ca varies more than any other bond length, and 
(ii) at least in simple cuprates with a single Cu02 layer 
in the unit cell it correlates with the maximal Tc - the 
longer ca, the higher Tc. At least, the first fact can 
be understood: ca is 'soft' because apical oxygen has no 
hard contact with the nearest copper ion; rather, it "levi- 
tates" on the electrostatic potential - a structural feature 
peculiar to certain layered oxides with alternating ionic 
planes of opposite charge [TT]. This makes apical oxygen 
prone to very large displacements - e.g., in HgBa2Cu06 
one finds ca ~ 2.8 A, longer by 0.9 A than the in-plane 
Cu-0 bond; coincidentally, this compound has the high- 
est Tc = 97 K among all singlc-Cu02-layer cuprates [l8] . 

It is thus important to find out what happens to the 
apical oxygen in LSCO-LCO bilayers; however, standard 
X-ray diffraction (XRD) is not well suited for this — one 
needs to "get inside the unit cell" and look for individ- 
ual atomic displacements. For this, the most suitable 
technique is the Coherent Bragg Rod Analysis (COBRA) 
method [HISnilSIJIllllllMj- However, COBRA is most 



effective for films that are just a few unit cells (UCs) 
thick, and in the case of HTS compounds fabrication of 
ultrathin films with bulk properties has been proven to 
be extremely challenging. Fortunately, we had a tech- 
nical solution at hand — a unique atomic layer-by-layer 
molecular beam epitaxy (ALL-MBE) system with proven 
capability of fabricating ultrathin HTS layers [25] |26] . 

For this study we have synthesized by ALL-MBE a 
number of (n x LSCO -I- m x LCO) bilayers, where (n,m) 
determine the thickness of the respective layers expressed 
as the number of UCs. In this paper, we show the CO- 
BRA results for two of these, (2.5,2.5) and (2,3). The 
films were deposited at T = 650°C and p = 9 x 10~^ 
Torr of ozone and subsequently cooled down under high 
vacuum to drive out all the interstitial oxygen. We used 
10 X 10 X 1 mm ^ single-crystal LaSrA104 (LSAO) sub- 
strates polished with the large surface perpendicular to 
the (001) direction. The substrate lattice constants are 
ao = bo = 3.755 A, Co — 12.56 A; the films are pseudo- 
morphic with LSAO and under compressive strain. The 
crystal structure of LCO is illustrated in Fig. la. Atomic 
force microscopy scans over a large (10 x 10 fini^) area 
showed root-mean-square surface roughness of 0.25 nm 
in the (2,3) and 0.11 nm in the (2.5,2.5) bilayer sam- 
ple; this is significantly less than the 1 UC step height 
which in LSCO is 1.32 nm. Magnetic susceptibility was 
measured via two-coil mutual inductance technique and 
revealed sharp superconducting transitions at Tc = 34 K 
in the (2,3), and Tc = 36 K in the (2.5,2.5) bilayer, signif- 
icantly higher than the values reported for (n,m) bilayers 
in Rcf . [3] , which is remarkable given that these films are 
only 5 UC thick. This was also confirmed by measuring 
the electric resistance (see Fig. lb) after the X-ray scat- 
tering experiments were completed and gold pads were 
evaporated to enable four-point-contact measurements. 

The atomic structure of the LSCO/LCO bilayer film 
was investigated at beamline ID-33 of the Advanced Pho- 
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FIG. 1: (Color online) A simplified structure model (one- 
half the crystallographic unit cell) of La2Cu04 and the trans- 
port property for a bilayer film, (a) At room temperature, 
the structure is tetragonal and the space group is I4/mmm. 
Noted that the La(Sr)-0 layers are strongly corrugated, exag- 
gerated in this sketch for clarity, (b) The electric resistance 
of the (2.5,2.5) bilayer, measured by the four-point-contact 
technique, as a function of temperature. Inset: a schematic 
of the bilayer on a LSAO substrate. 



ton Source by measuring the diffraction intensities along 
the substrate-defined Bragg rods. The sample and a PI- 
LATUS 100k photon-counting pixel detector |27j were 
mounted on a six-circle goniometer in Kappa geometry. 
The experimental set-up and procedures were described 
in detail in previous works [SHI HI] . Ten symmetry in- 
equivalent Bragg rods were recorded with a maximum 
value for the vertical reciprocal space coordinate of Lmax 
= 10.5 r.l.u. (reciprocal lattice units) and a sampling 
density of 50 points per r.l.u.. The X-ray flux after the Si 
(1,1,1) monochromator crystal was 3x10^^ photons/sec 
at a wavelength of A = 0.8266 A. For all Bragg rod mea- 
surements, except for the (0,0, L) rod, the angle of inci- 
dence had a fixed value of 3.5°. The X-ray beam was 
focused to 0.1 mm (V) x 0.2 mm (H), resulting in a 2 
mm long X-ray footprint. The background and diffuse X- 
ray scattering contribution were removed efficiently and 
accurately using the PILATUS detector images. The fi- 
nal results were then normalized by taking into account 
the beam polarization and Lorentz factors. The results 
were subsequently analyzed using the COBRA method 
[H [511 [21 1231 [^. In general, COBRA uses the mea- 
sured diffraction intensities and the fact that the com- 
plex structure factors (CSFs) vary continuously along the 
substrate-defined Bragg rods to determine the diffrac- 
tion phases and the CSFs. The CSFs are then Fourier 
transformed into real space to obtain the 3-dimensional 
electron density of the film and the substrate with sub- 
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FIG. 2: (Color online) Representative Bragg rods of 
LSCO/LCO system (open diamond) and calculated diffrac- 
tion intensity obtained from COBRA-determined electron 
density (solid line). 

Angstrom resolution. 

The experimental data of three representative Bragg 
rods are shown in Fig. 2(a). Notice that the diffrac- 
tion intensity along the Bragg rods (excluding the Bragg 
peaks) vary over more than 4 orders of magnitude with 
excellent signal-to-noise ratio. The reference structures 
chosen as the starting point for the COBRA analysis were 
the bulk LSAO structure and the tetragonal LSCO /LCO 
bilayer with the nominal layer atomic positions. In our 
numerical simulations, the topmost 4 UCs of the sub- 
strate were allowed to deform, however the resulting 
deformations turned out to be very small. The CO- 
BRA method uses the approximation that at two ad- 
jacent points along the Bragg rod the change in CSFs 
contributed by the unknown part of electron density is 
negligible compared to the change in CSFs contributed 
from the reference structure [20 j . The use of this ap- 
proximation allows COBRA to converge very quickly to 
approximately the right solution but not to the exact one. 
To overcome this limitation wc further refined the CSFs 
using the Difference-Map algorithm introduced by Elser 
P5] and recently applied to thin films [53]. Using the 
COBRA solution as the starting point for the Difference- 
Map algorithm and using a proper filter program that 
takes advantage of the fact that the CSFs vary continu- 
ously along the Bragg rods, the Difference-Map algorithm 
converges after about 20 iterations; the convergence ac- 
celerates by about two orders of magnitude. As seen in 
Fig. 2, the final calculated and measured intensities are 
in very good agreement. Similar agreement was found 
for all other Bragg rods and the overall X-ray reliability 
factor R = ^" = 0.02; here, Fq and F^ are the ob- 

served and the calculated diffraction amplitudes, respec- 
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FIG. 3: (Color online) The electron density variation, de- 
termined by COBRA, along the [0, 0, Z] column of atoms 
as indicated by the dashed line in Figure la. Note that the 
topmost four unit cells of the substrate are included in the 
structure refinement. The left and right dashed lines repre- 
sent the nominal LSAO/LSCO and LSCO/LCO interfaces, 
respectively. 



tively. To the best of our knowledge, so far there has been 
only one attempt to determine the structure of a thin film 
using the Difference Map method [29]. In that study, the 
atomicity constraint was imposed and over 2,000 iter- 
ations were needed to achieve convergence. Our anal- 
ysis shows that the combined COBRA/Difference Map 
method combines the best features of both methods and 
ensures rapid convergence to the correct solution without 
the need to use the atomicity constraint. 

The CSFs obtained have been Fourier transformed into 
real space yielding the 3D electron density (ED). As an 
example, we show in Fig. 3 the ED of a (2,3) bilayer sam- 
ple along the [0,0, Z] line that goes through the La(Sr), 
O, and Cu(Al) atoms. Two points should be stressed. 
First, as seen the ED has almost no negative parts. To- 
gether with the excellent agreement between the calcu- 
lated and measured diffraction intensities, this suggests 
that the ED is very close to the correct one. Second, all 
the atoms including the oxygens can be clearly identified 
and their positions determined with sub- Angstrom reso- 
lution. The small ED intensity fluctuation below -53 A 
provides a measure of the inaccuracy in the ED and as 
seen it is small even compared to the oxygen ED. 

The atomic positions in the Z direction were accu- 
rately determined by fitting a Gaussian to each peak. 
We determined the size of the unit cell in the Z direction 
by measuring the distance between consecutive pairs of 
La(Sr) and Cu atoms. The results are shown in Fig. 4 
inset. Each point corresponds to an average of 4 La-La 
distances and 2 Cu-Cu distances. The measured lattice 
constant of the bilayer film is 13.304 ± 0.016 A and is 
larger by 0.148 A than that of the bulk LCO (cq = 13.156 

A). 

While the changes observed in the unit cell sizes are as 



3.0 
2.8 
2.6 
2.4 
2.2 
2 
1.8 



(a) 



♦ (2,3) bilayer 

• (2.5,2.5) bilayer 
A Average 



-4 -3 -2 -1 1 2 3 4 5 
Z position (unit cells) 



i 




-4 -3 -2 -1 1 2 3 4 5 
Z position (unit cells) 



FIG. 4: (Color online) Evolution of the inter-atomic distances, 
(a) The measured Cu(Al)— apical O distance, ca, varies as a 
function of the nominal position of Cu(Al) atoms inside the 
refined structure. The data from two representative bilayer 
samples and the average over the two are presented. The lower 
and upper arrows represent the bulk values of ca for LSAO 
and for LCO, respectively, (b) The comparison of ca, ci, and 
C2, averaged for each unit cell, as a function of Z position. 
Inset: the lattice constant co as a function of Z. The dotted 
line represents the bulk LSAO value. The horizontal dashed 
line is the average value of Co in bilayers extracted from the 
electron density, as described in the text. In both (a), (b) 
and the inset, the vertical dashed lines represent the nominal 
LSAO/LSCO interfaces, respectively. 



expected from the strain and the elastic parameters |30j , 
the variations in Cu-apical O and La-apical O distances 
are quite unexpected. The distances Cyi,Ci and C2 are de- 
fined in Fig. la; the distances labeled c^,Ci and C2 would 
be their symmetry equivalents in bulk samples, but in 
thin films they could differ in principle. For our samples, 
the measured values for the primed and unprimed dis- 
tances were in fact equal within the experimental error, 
except at the LSAO/LSCO interface. The measured val- 
ues averaged over ca and c^ are shown in Fig. 4a. The 
diamond and circular dots represent the distances mea- 
sured in (2,3) and (2.5,2.5) bilayers, respectively. The 
triangular dots are averages over the two samples. Ev- 
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ery pair of triangular dots corresponds to one UC. The 
dashed vertical lines represent the nominal LSAO /LSCO 
interfaces. The arrows on the right indicate ca as mea- 
sured in the bulk samples. The results show that, within 
the experimental error, the values of ca in the substrate 
are equal to those in the bulk but they are very different 
in the film. In both LSAO and LCO bulk crystals [TZ] . 
Ca is equal to 2.41 A. In the metal layer closest to the 
substrate Ca — 2.3 A, and it then rises steadily all the 
way to Q-A = 2.75 A — a change of 0.45 A. 

In Fig. 4b we display ca as well as the La-apical O 
distance, C2, and the La-Cu02 plane distance, Ci. Each 
point represents an average over the two bulk-symmetry- 
equivalent distances and over the two measured samples. 
Notice that Ci changes by less than 0.1 A. On the other 
hand, ca increases by about 0.45 A, while C2 decreases by 
about 0.25 A. Close to the film surface, the apical oxygen 
atoms are displaced away from the nearest Cu atoms. 
The La atoms are displaced towards the closest Cu02 
plane, but by a smaller amount, while the separation 
between two adjacent Cu02 planes remains constant. 

According to Ref. 15, ca — 2.7 A should correspond to 
a Tj, of 80 K at the optimum doping. However, from Ref. 
31 we know that the hole density drops sharply on the 
I side of the interface and the screening length is equal 
to 6 ± 2 A. This implies that on the I side and next to 
the M-I interface only one or two Cu02 layers are doped 
via carrier accumulation while the others remain insulat- 
ing. Thus, unfortunately, we have a mismatch: in the 
optimally-doped LCO layer, ca is close to its standard 
(bulk) value, while it is greatly elongated only in insu- 
lating layers. It is tempting to speculate that one could 
create LSCO-based samples with Tj, much higher than 36 
K, perhaps as high as 80-90 K, if only one could achieve 
Ca elongation and optimal doping in the same LCO layer. 
An obvious avenue for further research is to try making 
I layers even thinner, thus bringing the interface super- 
conductivity closer to the film surface. Another is to 
try engineering more sophisticated hetero-structures and 
superlattices combining LCO with other metallic oxides 
(nickelates, zincates, etc.). 

In summary, we have used ALL-MBE to synthe- 
size precise ultrathin bilayers using metallic but non- 
superconducting LSCO and insulating LCO blocks, 
and observed interface superconductivity with Tc = 
34-36 K, significantly higher than before. We have 
used synchrotron X-ray diffraction and the combined 
COBRA/Difference-Map phase-retrieval method to de- 
termine accurately the atomic structure and found the 
unit cell size to be constant despite dramatic atomic dis- 
placements within the cell. In particular, near the surface 
the Cu— apical O increases greatly, by as much as 0.45 A, 
while it is known that variations in the apical oxygen po- 
sition strongly affect T^. We conclude that in cuprates 
the crystal structure can be modified in near-surface lay- 
ers, and in such a way that superconductivity properties 



can be dramatically altered. This result amplifies the im- 
portance of high quality surface structure determination 
in conjunction with surface sensitive probes of electronic 
states such as scanning tunneling microscopy or angle- 
resolved photoemission spectroscopy. 
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